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Foreword 


This final report is submitted by SECA, Inc. , to the National 
Aeronautics and Space Administration, George C. Marshall Space 
Flight Center, Alabama, for Contract NAS8-38243, entitled "Orbital 
Maneuvering Vehicle (OMV) Plume and Plume Effect Study . This 
report presents the results of work performed by SECA to fulfill 
the requirements of Contract NAS8-38243. 


Two other documents already published which are summarized in 
this report also comprise part of the final report. These documents 
which are contained as appendices to this report are: 


"Definition of the MRE15 Engine Exhaust Plume Flowfield", 
SECA-TR-90-03 , 30 March 1990 

"Definition of the OMV GN2 Thruster Exhaust Plume Flowfield", 
SECA TR-90-13 , 15 August 1990. 
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Summary 

The object of this contract was to characterize the Orbital 
Maneuvering Vehicle (OMV) propulsion and attitude control system 
engine exhaust plumes and predict the resultant plume impingement 
P^’®^®vre, heat loads, forces and moments. The resultant 
environments could then be used by the designer to insure adequacy 
of the structural, attitude control and thermal control systems. 

The original contract provided for two increments of funding 
for the entire study. However, late in the first year of this 
study the OMV project was terminated so that the second increment 
was not funded. As a result of the reduced funding, not all tasks 
were completed. The following studies were completed: 

1. The MRE15 (main propulsion and primary attitude control 

axisymmetric exhaust plume was characterized. 

2. The Gaseous Nitrogen (GN2) secondary attitude control 

axisymmetric exhaust plume was characterized. 

3. An axisymmetric Monte Carlo simulation of the MRE15 
exhaust plume was performed. 

4. Initial estimates of the base heating environments to the 
propulsion model were determined for an axisymmetric 
Monte Carlo MRE15 plume. 

Plirop code modifications and input files for locating 
the 28 hydrazine attitude control thrusters, 24 GN2 
attitude control thrusters, and four propulsion module 
thrusters were generated for input to the plume 
impingement program. These files could be used to 
perform the impingement environment prediction for the 
OMV. 

, The techniques, methods and models for predicting plume 
induced environments presented in this report can be applied to the 
Cargo Transfer Vehicle (CTV) which is presently being considered 
as part of the NLS vehicle. 
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1 . 0 Introduction 

This final report describes the results of work performed by 
personnel of SECA, Inc. for the NASA/MSFC Induced Environments 
Branch (ED33) under Contract NAS8-38243. This incrementally funded 
contract entitled, "Orbital Maneuvering Vehicle (OMV) Plume and 
Plume Effects Study", contained three major tasks which were to be 
completed if the incrementally funded contract were fully funded. 

Task 1 of the contract called for analyzing the Propulsion 
Model (PM) engine combustion chamber, nozzle and exhaust plume 
flowfield using continuous methodologies. This effort is 
described in Section 2.1 of this report. 

Task 2 of this contract called for the analysis of the PM 
engine flowfields using axisymmetric and three dimensional 
Monte Carlo techniques. The axisymmetric Monte Carlo analysis 
was performed and is described in Section 2.2 

Task 3 of the contract called for the analysis of the cold gas 
gaseous nitrogen (GN2) attitude control thruster, hydrazine 
attitude control thruster, exhaust plumes and the plume 
induced environments due to the attitude control thrusters. 
The plume analysis for this task was completed and part of the 
impingement analysis was performed. These results are 
presented in Section 2 . 3 
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2.0 Technical Discussion 

2.1 Propulsion Module Axisymmetric Plume Model Development 

The OMV base line design prior to and at the beginning of this 
contract utilized four Variable Thrust Bipropellant Engines (VTE) 
as the primary orbit change propulsion system on the Propulsion 
Module (PM) . As the design matured the VTE engines were replaced 
with MRE15, 15 lbf Hydrazine Thrusters. These engines are the same 
engines which are used as part of the attitude control thrusters. 
The continuum analysis of this engine is discussed in Section 2.3 
of this report. The continuum analysis of the PM engine/plume 
flowfield was necessary to provide the input data to perform the 
axisymmetric Monte Carlo analysis of the PM flowfield in order to 
calculate the base heating environment to the OMV. The 
axisymmetric Monte Carlo analysis and results are described in 
Section 2 . 2 

2.2 Axisymmetric PM Plume Direct Simulation Monte Carlo (DSMC) 

Analysis 

The object of this task was to provide an estimate of the OMV 
PM base environment due to a single engine plume. Figure 1 shows 
the PM engine layout for the OMV. There are four MRE15 engines 
located in the base as shown in Figure 1. When these engines fire 
the plumes interact and cause flow reversal into the base which 
results in heat and pressure loads which need to be accounted for. 
Due to the relatively low density of the flow the only way (lacking 
applicable experimental data) to predict the interaction of the 
flowfields is to perform a three dimensional analysis using Direct 
Simulation Monte Carlo (DSMC) techniques. The approach for 
performing this analysis required three steps. Step 1 requires 
performing a continuum analysis of the MRE15 engine and plume. 


2 



SECA-TR-92-02 


This analysis is described in Section 2.3. Step 2 uses the 
continuum analysis as input and performs an axisymroetric DSMC 
analysis. Finally, Step 3 uses the results of Step 2 to perform 
a full three dimensional Monte Carlo analysis to determine the base 
environment. Since the second increment was not funded, only the 
axisymmetric analysis was performed. 

DSMC methods have become more and more popular for solving 
complex flowfields in low density flows where transitional non- 
equilibrium and specie diffusion are important. Three dimensional 
DSMC methods could be used to solve the entire 4 engine plume 
flowfield starting at the exit plane of the engine but the computer 
time necessary to do this would be excessive. To reduce the 
computer resource requirements an axisymmetric DSMC calculation of 
a single MRE15 engine was performed An axisymmetric DSMC MRE15 
calculation was used to characterize the flowfield in the near 
region around the engine including the engine heat shield. The 
axisymmetric results were then used to calculate the heat loads to 
the engine heat shield region. 

The basic code which was used to perform the axisymmetric DSMC 
MRE calculations was the GAPS code (Ref. 1). This model was 

developed to solve unique design problems for the Galileo 
Spacecraft. 

The direct simulation Monte Carlo method was originated by G. 
A. Bird and has been developed by numerous modelers over the past 
decade. This method is a computational technique for the modeling 
of a real gas flow at a molecular level with several thousand 
simulated molecules. The velocity components and position 
coordinates of the simulated molecules are stored in the computer 
these are modified with time as the molecules are followed 
through representative collisions and boundary interactions in the 
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simulated flowfield. 

The initial state of the gas is specified and boundary 
conditions corresponding to the desired flow are applied at time 
zero. The flow is always unsteady, but the boundary conditions 
are generally such that a steady flow is obtained as the large-time 
state of the unsteady flow. Time is advanced by discrete steps of 
magnitude 10 At which is small compared with the mean free time 
between collisions. The simulated flowfield is divided into a 
network of spatial cells of such size that the change in flow 
properties across each cell is small. All macroscopic quantities 
of interest are obtained as averages of the randomly-fluctuating 
molecular states within the various cells. 

The molecular motion and the collision processes are uncoupled 
over the time interval At by the repeated application of the 
following procedure: 

(i) All the molecules are moved through distances appropriate 
to their instantaneous velocity components and At. 

(ii) A representative set of collisions, appropriate to At, 
is computed among the molecules, and the pre-collision 
velocity components of those molecules involved in the 
collisions are replaced by their post-collision values. 

Therefore, along any molecular path, collisions can only occur at 
intervals that are integral multiples of At. Since At is small 
compared with the mean free time, the resulting distortion of the 
path is small. 

The direct simulation procedure is equivalent to a numerical 
solution of the Boltzmann equation. The left hand side of that 
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equation is represented directly by the process of following the 
trajectories of the individual molecules. The collision integral 
on the right hand side is properly represented by sampling 
collision partners from the collision probability appropriate to 
the particular molecular model, and by calculating collisions at 
a rate consistent with the local collision frequency. It can be 
shown that, as the uncoupling interval At and the cell size are 
made very small, and the number of molecules is made large, the 
Monte Carlo simulation differs from the Boltzmann equation only by 
random statistical fluctuations. 

Since the conservation equations of continuum fluid mechanics 
are embodied in the Boltzmann equation, solution of the Boltzmann 
equation gives the proper description of gas flows at all density 
levels. Therefore, the direct simulation technique can, in 
principle, be used to describe the entire rarefaction spectrum from 
the collisionless to the continuum limit. The practical 
constraints for the application of this method to continuum flows 
are the limited computer resources; memory and time. For this 
reason, the GAPS program has been designed to only model the flow 
outside the "core" of the exhaust plume. The results of the 
continuum MRE15 analysis serve as the specified boundary for the 
GAPS code. Aerospace Corporation has performed DSMC studies (Ref. 
2) using the basic GAPS code. Earlier studies (Ref. 3) identified 
some areas of improvement in the original Bird (Ref. 4) formulation 
which have been included in the version of the GAPS code which was 
used for this study. These differences are confined to inflow 
boundary conditions, spatial cell construction and the referencing 
of molecules to a particular cell. 

Inflow boundaries must generally be located within a flow 
where continuum solutions are still valid. It is then possible to 
derive an explicit relation for the one-way flux of molecules 
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through an arbitrarily oriented plane element within this flow. 
The net inflow is the difference between the inward and outward 
one-way flows and can be expressed as 

N n = N + - N' 

Bird’s approach is to calculate N + and to introduce an essentially 
fixed number of molecules per time step through the boundary 
element (fluctuations occur because only an integral number of 
molecules can be introduced in any one time increment) . The 
outgoing flux at any element is then established by the collision 
processes within the cells adjacent to the boundary and this in 
turn also establishes the net flow. It was found that this 
approach did not reproduce the average net flow to the accuracy 
desired in the simulation of nozzle and plume flowfields. To 
overcome this a fixed net flow boundary condition is used and the 
number of inflow molecules is determined from 

N + = N n + N" 

This requires a small amount of additional storage because the 
value of N must be determined for every inflow element during a 
given time step. The procedure ensures that the various molecular 
fluxes at the boundary are matched to the continuum results. 

The spatial structure of cells is based upon the construction 
of distinct regions each of which is further subdivided into 
computation cells. Each region allows independently selected time 
steps and molecular weighing factors. These features are crucial 
to simulations which involve as much as, or more than, ten order 
of magnitude variations in basic flow properties such as density 
or in controllable features such as cell dimensions. 

The basic geometry of a region is a closed quadrilateral which 
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region. Then the results of the lip calculation were supplied to 
the complete axisymmetric DSHC model as input conditions much th 
same as the continuum plume source segments. By calculating 
lip region alone, then patching it into the complete solution, 

considerable computer resources were saved. 

The GAPS code has been used to model the Voyager spacecraft 
plumes. This model has been compared against flight calorimeter 

data with very good results. 

Three DSMC cases were run for the single MRE15 axisymmetric 
plume, case 1 was a free plume without an OMV base heat shield. 
Case 2 had a heat shield .375 inches ahead of the lip while case 
3 had the heat shield 1 inch ahead of the heat shield. Cases 2 and 
3 were calculated to provide a sensitivity of peak heat rate an 
distribution to the distance the engine projected out from e 

base. 

Figure 3-6 present the DSMC results for the MRE15 case 1 that 
has no vehicle base. Figure 3 presents the velocity vectors 
Figure 4 the density contours, Figure 5 rotational temperature 
contours and Figure 6 the translational temperatures It is 
apparent that the plume goes out of continuum at relatively low 
angles off the starting surface by examining the rota lona 
translational temperatures. 

Figures 7 and 8 present the results of case 2 where the base 
is .375 inches upstream of the nozzle exit. Figure 7 s ows 
velocity vectors, and Figure 8 shows density contours. Figure 
shows density contours for case 3 where the base is 1 inch upstream 

of the lip. 

AS was previously mentioned, the purpose of running cases 2 
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and 3 was to get an estimate of the base heating environments for 
a single engine as a function of nozzle exit plane location with 
respect to the base. Figure 10 presents the heating rate 
distribution to the base for cases 2 and 3 as a function of 
distance from the nozzle centerline. By examining Figure 10, it 
is obvious that the closer the exit plane is to the base the higher 
the heating rate. The peak heating rate for the .375 inch spacing 
is .029 BTU/Ft 2 /sec while the peak for the 1 inch separation plane 
is .01 BTU/Ft 2 /sec . The location of the peak heating rates are at 
different positions on the base although in general fairly near the 
lip. These results are qualitatively the same as what would result 
from a continuum impingement analysis, although the continuum 
results would be two to three times higher. 

A three dimensional analysis of the four engine thruster would 
produce higher peak heating rates than the single engine due to 
plume interaction and the reversal of more mass towards the base. 
An estimate of the peak heating on the order of . 1 to . 3 
BTU/Ft z /sec depending on the engine spacing and location of the 
exit plane. 

The results of cases 1-3 and the program used to perform 
these calculations are resident on the MSFC EADS computer system. 
They can be obtained by contacting Y. C. Lee at ED33. 

The three input files are contained on the first file of the 
magnetic tape that accompanies this final report (See Section 2.4) . 
The files are: CASE1.INP, CASE2 . INP and CASE3.INP. The magnetic 
tape was generated in a VAX/VMS Backup format at 6250 BPI . 

2.3 Hydrazine and Cold Gas Thruster Plume and Plume Impingement 

Modeling for the OMV Vehicle. 
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This section describes the modeling of the MRE15 and cold gas 
nitrogen attitude control plume calculations. Additionally, the 
input data files that were generated for the plume impingement 
induced environments are described. Not all plume impingement 
input data were generated since the second contract increment was 
not funded. None of the plume impingement environments were 
generated. The plume impingement program (PLIMP-Ref. 5) which was 
to be used for performing this analysis is also described along 
with the modifications which were made to it to perform the OMV 
Plume Impingement Induced Environments. 

2.3.1 MRE15 Hydrazine Plume 

The MRE 15 is a 15 lbf hydrazine engine which was to be used 
as the primary engine for the Reaction Control System (RCS) as 
well as the thrusters for the Propulsion Module (PM) of the 
Orbital Maneuvering Vehicle (OMV) . The MRE 15 plume character- 
ization is required in order to predict the MRE 15 plume induced 
environments to the OMV itself, as well as to payloads it is 
retrieving or deploying. 

The analysis of the MRE 15 nozzle and plume flowfields was 
performed using the Reacting and Multi-phase Computer Code RAMP2 
(Ref. 9) . RAMP2 is a method of characteristics based code used 
to solve inviscid supersonic flowfields. 

The purpose of the RAMP2 code is to solve a rocket exhaust 
flowfield in sufficient detail so that the plume may be used to 
adequately specify design environments for systems which are 
influenced by the plume. For high altitude plumes, an adequate 
definition of the plume requires the inclusion of chemistry and 
nozzle boundary layer effects. The TRAN72 (Ref. 7) program and 
BLIMPJ (Ref. 8) codes are included in the RAMP2 code along with 
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the RAMP2F flowfield code to treat chemistry and boundary layer 
effects . 

In order to use the original RAMP2 code (Ref. 9) to solve a 
high altitude plume, the following steps are required. First, 
the TRAN72 program input data is prepared and executed to 
generate a data tape describing the thermodynamic characteristics 
of the post-combustion gases. Next the RAMP2F flowfield data are 
prepared and the nozzle flowfield is solved using the TRAN72 
program thermodynamic data tape as input. Then in order to 
adequately describe the nozzle boundary layer, the BLIMPJ code is 
executed using an input data file and flowfield tape generated by 
the RAMP2F nozzle solution. Finally, the exhaust plume is 
generated by using the nozzle solution and boundary layer 
solution to generate an exit plane start line that is used to 
initiate the plume solution. Thus, the generation of a high 
altitude plume could require up to four different executions of 
programs (TRAN72 , RAMP2F, BLIMPJ, and RAMP2F) for the 
specification of the most detailed and accurate results. 

Physical input data are required only for the TRAN72 and first 
RAMP2F execution. All data required for the BLIMPJ code and 
second RAMP2F execution are generated internally to the program 
and/or communicated via temporary files. 

Recently RAMP2 has been modified (Ref. 6) which will better 
handle nozzles which operate at a low chamber pressure and/or 
high area ratio. Motors of these types may have boundary layers 
which contain a significant percentage of the mass flow of the 
motor. Large boundary layers can interact with the inviscid 
portion of the flow and alter the flowfield from what would be 
calculated, assuming a pure inviscid calculation. The original 
version of RAMP2 which was mentioned above (Ref. 9) could not 
account for the interaction of the boundary layer and the 
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inviscid flowfield. RAMP2 made a single pass through the nozzle 
and performed a single boundary layer solution using the BLIMPJ 
code. These solutions were then merged at the exit plane and 
plume solution was performed. 

The latest version of RAMP2 (Ref. 6) was used to calculate 
the MRE 15 plume. This version allows a first order 
approximation of the interaction of the boundary layer and 
inviscid flowfield. The basic method involves making two passes 
through RAMP2 to calculate the nozzle flowfield and two passes 
through the BLIMPJ code to calculate the boundary layer. The 
first pass through the RAMP2 code calculates the nozzle flowfield 
using the original nozzle contour. After the first nozzle 
solution, the BLIMPJ code is executed using the conditions which 
exist at the nozzle wall (from the first nozzle solution) as the 
inviscid streamline boundary layer edge conditions. The results 
of this calculation including the boundary layer displacement 
thickness as a function of nozzle station are stored on FORTRAN 
Unit 2. RAMP2 is then executed again and the nozzle contour is 
modified using the displacement thickness information contained 
on Unit 2. After the modified contour flowfield solution is 
completed, BLIMPJ is again executed to calculate the nozzle 
boundary layer using edge conditions taken at the wall of the 
modified nozzle flowfield solution. While in the BLIMPJ 
solution, a check is made once the boundary layer calculations 
have been completed. If the pressure at the lip of the nozzle 
varies by more than 10% from the pressure at the inside edge of 
the boundary layer in the inviscid solution (modified contour) 
then another pass is made through BLIMPJ with the edge conditions 
interpolated from the nozzle solution based on the location of 
the boundary layer edge in the flowfield. These two passes are 
made within the last BLIMPJ execution and require no intervention 
by the user. This second pass is made to ensure that the 


12 



SECA-TR-92-02 


boundary layer and inviscid flowfield properties match at the 
edge of the boundary layer. The RAMP2 , TRAN72 and BLIMPJ 
versions which were used for this analysis are contained in File 
2, 3 and 4 of the accompanying magnetic tape. 

Recently, a model has been developed which fits the RAMP2 
flowfield results using a source flow approximation. This model 
called SFPGEN (Ref. 10) was developed to support a Space Station 
Contamination Model developed by SEA, Inc. for NASA (Ref. 11 ) . 

The source flow plume model generator is a code which 
determines the flow properties from a RAMP2 flowfield and 
generates a formatted file that can be used by the PLIMP code to 
rapidly determine plume impingement pressures and heating rates 
out to large distances from the exit plane. This model was used 
to calculate the MRE15 source flow plume. A more detailed 
description of the source flow model can be found in Appendix A. 


The source flow code and input files for the MRE15 source 
flow plume are on file 5 of the magnetic tape. The complete 
description of the MRE15 analysis and results can be found in 
Appendix A (Ref. 12). File 6 of the magnetic tape contains all 
MRE15 input/output files for the MRE15 RAMP2 analysis as well at 
the source flow plume file (MRE15UN16.DAT) 

2.3.2 OMV GN2 Plume Analysis 

The OMV GN2 RAMP2 flowfield analysis and source flow plume 
analysis were performed using the same tools that were described 
in Section 2.3.1. The only difference in the analysis was that 
only a single pass was made for the boundary layer analysis and 
ideal gas was assumed. The OMV GN2 nozzle has an area ratio of 
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elemental area local flow regime is determined by comparing the 
mean free molecular path of the local plume to the projected body 
dimension. In the continuum flow regime, the local impact 
pressures are computed using modified Newtonian theory which 
includes a correction factor to account for the higher pressures 
noted experimentally for small incidence angles (Ref. 16) . In 
the free molecular flow regime, the body local flow properties 
are assumed to depend only on collisions of the gas molecules 
with the body surface. Environmental equilibrium is assumed so 
that a Maxwell ian-type flow exists and the body local flow 
properties are calculated from kinetic theory and Maxwell's law 
of distribution of molecular velocities. The local force is then 
computed from the sum of the forces produced by the incident and 
reflected molecules. In the transitional flow regime the local 
impact pressure is calculated empirically as a function of the 
local Mach number and Knudsen number to yield a smooth transition 
from the continuum to the free molecular flow regime. 

This program contains several methods for calculating 
convective heat transfer rates to bodies immersed in rocket 
exhaust plumes. Each method is applicable to a particular 
situation, determined either by geometry or by the assumption 
made in obtaining the local flow properties. No single method 
currently will suffice for all shapes and flow conditions. In 
fact, a different method must be frequently used for the various 
subshapes of an impinged body. The method to be used on a given 
subshape must be specified in the input data. 

Several options are also available to the user for obtaining 
the local flow properties for use in the heat transfer analysis. 
In all of the methods, the local pressure depends on the angle 
that the local surface makes with the free stream flow. 
Fortunately, this method is applicable to a majority of 
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hypersonic flow problems. In the leading edge and corner 
regions, corrections have to be made to the pressure and heating 
values as calculated by the program. Free molecular heating 
rates are also calculated at each elemental area from an energy 
balance type equation and the decision of which flow regime to 
consider is simplified by a printout of the applicable flow 
regime based on the local Kn. 

The program can currently analyze the impingement flow 
properties on six (6) basic subshapes. These are: 1) a surface 

described by a conical equation, i.e., cylinder, cone, etc.; 2) 
rectangular flat plate; 3) circular plate; 4) right triangle; 

5) airfoil; and 6) a body of revolution in which the surface 
curvature as a function of the longitudinal axis is described by 
a polynomial equation. (A maximum of fourth-order curve can be 
used. ) 

The original version of the code will only handle up to 10 
subshapes for a given engine impinging on a body. The OMV 
vehicle model requires using many more subshapes. While multiple 
computer runs could be made for each engine, to calculate the 
entire environment from a bookkeeping standpoint, it is much more 
convenient to perform the calculations in one computer run. The 
version of the PLIMP code which is being delivered with this 
report has been modified to treat 50 subshapes. 

For free molecular impingement the subshape modeling as far 
as how to break it down into elemental areas is not important 
since only the location and impingement angle relative to the 
engine (along with plume mass flux, temperature and surface 
temperature) are important. The present version of the PLIMP 
code utilizes a model which treats the impingement of molecules 
in the free molecular regime utilizing a diffuse or cosine 
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scattering assumption. A 
be found in Ref. 17 . The 
scattering. 


very good description of 
following sketch depicts 


the model can 
this type of 



in the case of diffuse scattering, the incident gas 
molecules are absorbed by the surface and are subsequently 
reemitted at angles which are independent of the angle of 
incidence. This diffuse or cosine scattering is most typical of 
y OW energy lases incident upon amorphous surfaces 

bH SUrf f 6S " ay be Clea " in spacecraft cleanliness sense, 
ut have adsorbed gases or surface oxidation present, in this 

case, the normal and tangential components of the incident 
molecular momentum are transferred to the surface, and some 
additional amount, dependent upon the surface temperature is 
imparted normal to the surface by the reemitted molecules. 

Maxwell first proposed that some fraction f of the incident 
molecules is absorbed by the surface and is subsequently 
reemitted diffusely, with a distribution of velocities 
corresponding to that in a still gas at the temperature of the 
surface The remaining fraction, ! - f, is perfectly reflected. 
All of the incident normal momentum and a fraction f of the 
incident tangential momentum are assumed to be transferred to the 
surface. Additional momentum is imparted normal to the surface 
by the diffuse reemission of the f fraction and the specular 
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reflection of the 1 - f molecules. Knudsen, in an attempt to 
characterize the energy exchange between gas and surface, 
introduced a coefficient of accommodation to describe the degree 
to which the incident molecules were "accommodated" to the 
temperature of the surface. Expressed in terms of energy flux 
per unit of surface area, this energy accommodation coefficient 

a e may be defined as 



where the subscript i refers to the incident energy flux, r 
denotes the energy flux carried away from the surface by the 
reflected molecules, and s refers to that flux which would be 
carried away from the surface if the reemitted gas molecules were 
in Maxwellian equilibrium at the temperature of the surface. 

Expressing Maxwell's scattering characterization in similar 
fashion leads to an expression for what has commonly become known 
as the tangential momentum accommodation coefficient a t 
(corresponding to f) 


- r , 


Oft = 


where ,, denotes that component of the incident momentum which is 
tangent to the surface, and is the tangential component of 
momentum imparted to the surface by the reemitted molecules;,, 
denotes the tangential momentum which is imparted to the surface 
by the reemission of molecules in Maxwellian equilibrium at the 
temperature of the surface, and is zero by symmetry. 
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Believing that the momentum exchange process is yet more 
complex. Bell and Schaaf (Ref. 18) introduced what may be termed 
a normal momentum accommodation coefficient, a n , defined by 
analogy to the first equation above as 



where P refers to the component of momentum imparted normal to 
the surface, and the subscripts are as defined above. 

At the present time the PLIMP code assumes all accommodation 
coefficients are unity. The models for force and heating rate 
presently in the PLIMP code are identical to those presented in 
Ref* 17, however there presently is no tangential force 
calculation in PLIMP. Depending on orientation and size of a 
surface, the tangential component can significantly effect the 
forces and torques that the OMV could experience due to free 
molecular impingement. The PLIMP code was modified to include 
the tangential component of force. Making the same assumption 
that the incoming particles have a Maxwellian distribution of 
molecular velocities about the macroscopic gas velocity, the 
incident tangential momentum component may be written: 

T X = -ypU 2 X COS{0}/(S7T 5 ) 

Here p and U are the incident gas density and velocity, 0 is the 
impingement angle. 

7 is the ratio of specific heats, and S is the molecular 
speed ratio, defined as the ratio of the macroscopic gas velocity 
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to the most probable molecular velocity: 

S = U/ ( 2RT) °' 5 

Other quantities are defined as: 

r) = S sin{6) 

X = exp {-r? 2 } + [1 + erf {»?)] 

Since all of the incident molecules are reemitted diffusely 
from the surface with a Maxwellian distribution of velocities 
corresponding to the temperature of the surface, the outgoing 
tangential component of momentum removed from the surface is 
zero. Thus only the incident tangential component was 
incorporated into PLIMP. This modification is included in the 
version of PLIMP which is contained on file 8 of the accompanying 
magnetic tape. 

2. 3. 3. 2 OMV Hydrazine/GN2 Thruster Location and PLIMP Input 
Files 


This section describes the OMV attitude control thruster 
geometry and inputs to the PLIMP code. The plume impingement 
calculations were not performed since the second contract 
increment was not funded. 

Figure 12 shows the OMV body fixed frame coordinate system 
which was used as the reference coordinate system for setting up 
the PLIMP code. Figure 13 shows a schematic of where the 
hydrazine thrusters are located on the OMV. The GN2 thrusters 
are not shown but are located adjacent to each hydrazine thruster 
except for hydrazine thruster A7 , B7 , C7 , D7 for which there is 
no corresponding GN2 thruster. Figure 14 shows a typical layout 
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Of the OMV orbital Replacement Unit (ORU, and the location of the 

the h d 1 Sh ° WS th£ l0Cati ° n and Action cosines of 

Table ^shor “ rUSterS . ln the PLIMP coordinate system while 
s e same information for the GN2 thrusters. 

The input files for all 28 hydrazine thrusters were 
generated for the PLIMP code. These data are 

genera^ 1 " T** ' *" p «* mj were 

e previously mentioned funding limitations. 

rrL 0t ° ata Contained O" Magnetic Tape for OMV Plume and 
Plume Impingement Analyses 

The programs used to perform this study as well as th. l * 

ana output files for generating the results discussed in this 
eport are contained on the accompanying magnetic tape. The tape 
was generated at 6250 BPI using a VAX/VMS format. Table 3 
summarizes the contents of the magnetic tape. 

3.0 Summary and Conclusions 

the o^ i en rePOrt reSUltE ° f 3 stud * to characterize 

the OMV engine exhaust plume flowfields and predict the resultant 

P ume induced pressure and heat loads to the vehicle structure 

Continuum methodology was used to calculate the MRE15 and GN2 ' 

the Cl5 P exT S 't T itl0naUy DSMC -i^ric calculations of 

the expected b l Pert °™ d a " d resultant estimates of 

e expected base heating to the Propulsion Model from a single 

engine were presented. Modifications to the Plume Impingement 
Program described as „ eU as ^ ^ ^ ^ p P J 

OMV attitude control thrusters a (■ano 

... . . ers< A ta P e accompanies this reDort 

were uTeTTnT ^ 'T**'** C ° deS and in P u Voutput files which 

used and generated during this study. All subtasks were not 
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completed due to program termination and the lack of the second 
funding increment. 

The codes and methodology delivered and described in this 
report are directly applicable to calculating plume induce 
environments for the National Launch System Cargo Trans er 

Vehicle. 
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Table 1 . 


OMV Hydrazine Thruster Locations and Directions in Body 
Fixed Frame 


Thruster 


Location (in) 


Direction Cosines 
(Thrust Vectors) 



X 

Y 

- 

X 

Y 

7 



c 4 ^ f i 

-63.44 

1.000 

0.000 

0.000 

A 1 
A2 

— 6 ij , 5 0 

-1.74 

-i 74 

54. 16 
63.39 

— 6 3 . j 9 
-54. 16 

0.000 • 
0.000 

-0.985 

-0.174 

0.174 

0.983 

A-.- 

A* 

* c 

As 

-3o. 74 
-4? . 33 
-49.33 

61.69 
68.3° 
54. 16 

-61.69 
-54. 16 
- .5 3 . v 9 

-0.707 

0.000 

0.000 

-0.500 

-0.174 

-0.985 

0 . 500 
0.985 
0. 174 

El 

b: 

z ■ 

u -V 

o 4 

u ■ 

E 3 
Bo 

- a 0 . 5 0 

- 1 . 7 4 
-3b. 74 
-47.33 

63. 44 
65.39 
54.16 
ci.6? 
54. 16 
65. 39 

59.20 

54.16 
65. ■-* 
61.69 
65.39 
5 4 . 1 a 

1 . 000 
0.000 
0.000 

-0.707 

0.000 

0.000 

0.000 

-0.174 

-0.9S5 

-0.500 

-0.985 

-0.174 

0.000 

-0.933 
-0. 174 
-0 . 500 
-0.174 
-0.953 

r 

c: 

-40.5 i 

_ i * 7 4 

- t 4 7 4 

-39.20 
-54. 18 
-65 . 3* 

63.4-4 

65.39 

34.16 

1 . 000 
0.000 
0.000 

0 . 000 

0.985 
0. 174 

0 . 000 

- 0 . 174 
-0.965 

C 4 

U Cr 

C* 

- 3 s • 7 4 
-49.33 

- 3 1 • 6 9 
-63.39 
-54.16 

ol.o? 
54, 16 

TO 

-0.707 

0.000 

0.000 

0.500 
0. 174 
0.935 

-0.500 
-0 . 9E5 
-0.174 

0 1 

-1.74 

-63.44 
-ob. 39 

-59.20 
-34. 16 

1.000 

0.000 

0.000 
0. 174 

0 . 0 0 0 

0.953 

D 2 
D3 
C 4 
r< c 

Do 

-1,74 

— -7 1 

- -w 0 • t * 

-49.35 

-49.35 

-34. 16 
-61. 69 
-34. it 
-63.39 

- i 8 . O'? 

- a 1 . 6 9 
-aS. 39 
-54.16 

0.000 
-0 . 707 

0 . 0 0 0 
0.000 

0.983 

0.500 

0.965 

0.174 

0 . 17 4 
0 . 5 0 0 

0. 174 
0.933 


, A - • 

63.44 

-59.20 

1.000 

0.000 

0.000 

A7 

B7 

C7 

D7 

— 0 * 'j ij ' 

-60 . 50 
-60. 5 0 
-60. 50 

59.20 

-63.44 

-59.20 

63. 4» 
59 . 20 
-63.44 

1.000 

1.000 

1.000 

0.000 

0.000 

0.000 

0. 000 
0 . 000 
0 . 000 
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TaMe 2- p"e d C °^ am G r Th — and Directions ln Body 


Thruster 


Location (in) 


Direction Cosines 


A 1 


A4 

A3 

H C 

E i 


tz 


Zb 


i_: 4 

C3 


U 

D4 



x Y 



i . ooo 

0 . 0 0 0 

0. 000 

o, ooo 

-o. 965 

0.174 

0 - U 0 0 

“0. 174 

0.9S5 

-0.707 

-0. 482 

0.517 

0.000 

-0. 174 

0.985 

0 . ooo 

“0.935 

0.174 

I . ooo 

0 . 0 0 0 

0.000 

0.000 

“0. 174 - 

0. 935 

U . 0 0 0 

-0.985 - 

0. 174 

“0.707 • 

-0.517 - 

0. 482 

0.000 ■ 

*0.985 - 

0.174 

0 .000 - 

■0. 174 

0.985 

i . ooo 

0. 000 o. ooo 

0 . 0 o 0 

0.985 -<: 

). 174 

0 . ooo 

0. 174 -o 

-.985 

“ 0 .707 

0.452 -0 

'.517 

0. ooo 

0.174 -0 

. 985 

0 . 0 0 0 i 

0.965 -0 

. 174 

1 • O 0 0 . 0 0 0 0 

. 000 

o. 000 c 

>.174 0 

. 963 

0 . 0 0 0 0 

'•965 0. 

, 174 

•0. 707 o 

'.517 0. 

4e: 

0.000 o 

.995 0. 

1 74 

0.000 o 

.174 0. 

985 
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Fig. 2 Plume Domains and Source Lines for DSMC Analysis 
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Fig. 3 Velocity Vector Distribution of MRE15 Flowfield Without a Base, 
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Fig. 4 Density Contours For MRE15 Flowfield Without a Base, 
Calculated Using DSMC Model 







Fig. 6 Translational Temperature Contours For MRE15 Engine 
Calculated Using DSMC Model 
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Fig. 9 Density Contours For DSMC MREl 5 Plume With the OMV Base, 
1.0 Inches Ahead of Lip 






3. 


Centroid of 
Grid Element 
(typical) 

Shuttle configuration used as example 
since representation of it requires use 
of many subshapes. 

Wing subshapes are NACA airfoil 0012-64. 
This can easily be changed to another air- 
foil geometry. 

Impingement data and heating rates obtained 
at centroid of each grid element. 



Fig. 11 - Sketch of Subshape System and Grid Schene Used in 

fbrce, tonent. Contamination and Heating Analysis 
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1.0 Forward 

This document presents the results of work performed by 
SECA, Inc. for NASA, Marshall Space Flight Center under Contract 
NAS8-38243 . The NASA contracting officer representative for this 
study is Mr. P. R. Sulyma, ED33. 

2.0 Introduction 

This report presents the flowfield characteristics of the 
MRE 15 engine exhaust plume. The MRE 15 is a 15 lbf hydrazine 
engine which will be used as the primary engine for the Reaction 
Control System (RCS) as well as the thrusters for the Propulsion 
Module (PM) of the Orbital Maneuvering Vehicle (OMV) . The MRE 15 
plume characterization is required in order to predict the MRE 15 
plume induced environments to the OMV itself, as well as to 
payloads it is retrieving or deploying. This report presents the 
MRE 15 flowfield characteristics and describes the methodology 
which was used to characterize the plume. 

3.0 Discussion 

This section describes the MRE 15 exhaust plume flowfield 
which was generated using the RAMP2 (Ref. 1) code and a new 
source flow model (Ref. 2) . 

3.1 RAMP2 MRE 15 Flowfield 

The analysis of the MRE 15 nozzle and plume flowfields was 
performed using the Reacting and Multi-phase Computer Code (RAMP2 
Ref. 1) . RAMP2 is a method of characteristics based code used to 
solve inviscid supersonic flowfields. 

The purpose of the RAMP2 code is to solve a rocket exhaust 
flowfield in sufficient detail so that the plume may be used to 
adequately specify some design environment for systems which are 
influenced by the plume. For high altitude plumes, an adequate 
definition of the plume requires the inclusion of chemistry and 
nozzle boundary layer effects. The TRAN72 (Ref. 3) program and 
BLIMPJ (Ref. 4) codes are included in the RAMP2 code along with 
the RAMP2F flowfield code to treat chemistry and boundary layer 
effects. 

In order to use the original RAMP2 code (Ref. 5) to solve a 
high altitude plume, the following steps are required. First, 
the TRAN72 program input data is prepared and executed to 
generate a data tape describing the thermodynamic characteristics 
of the post-combustion gases. Next the RAMP2F flowfield data are 
prepared and the nozzle flowfield is solved using the TRAN72 
program data tape as input. Then in order to adequately describe 
the nozzle boundary layer, the BLIMPJ code is executed using an 
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inpiit data file and flowfield tape generated by the RAMP2F nozzle 
solution. Finally, the exhaust plume is generated by using the 
nozzle solution and boundary layer solution to generate an exit 
plane start line that is used to initiate the plume solution. 

us, the generation of a high altitude plume could require up to 
four different executions of programs (TRAN72, RAMP2F, BLIMPJ 
and RAMP2F) for the specification of the most detailed and 
accurate results. Physical input data are required only for the 
TRAN72 and first RAMP2F execution. All data required for the 
BLIMPJ code and second RAMP2F execution are generated internally 
to the program and/or communicated via temporary files. 

Recently RAMP2 has been modified (Ref. 1) which will better 
handle nozzles which operate at a low chamber pressure and/or 
area ratio nozzles. Motors of these types may have 
boundary layers which contain a significant percentage of the 
mass flow of the motor. Large boundary layers can interact with 
the mviscid portion of the flow and alter the flowfield from 
what would be calculated, assuming a pure inviscid calculation. 
The original version of RAMP2 which was mentioned above (Ref. 5) 
could not account for the interaction of the boundary layer and 
the inviscid flowfield. RAMP2 made a single pass through the 
nozzle and performed a single boundary layer solution using the 
BLIMPJ code. These solutions were then merged at the exit plane 
and plume solution was performed. 

The latest version of RAMP2 (Ref. 1) was used to calculate 
the MRE 15 plume. This version allows a first order 
approximation of the interaction of the boundary layer and 
inviscid flowfield. The basic method involves making two passes 
through RAMP2 to calculate the nozzle flowfield and two passes 
through the BLIMPJ code to calculate the boundary layer. The 

P ass through the RAMP2 code calculates the nozzle flowfield 
using the original nozzle contour. After the first nozzle 
solution, the BLIMPJ code is executed using the conditions which 
exist at the nozzle wall (from the first nozzle solution) as the 
inviscid streamline boundary layer edge conditions. The results 
°f this calculation including the boundary layer displacement 
thickness as a function of nozzle station are stored on FORTRAN 
Unit 2. RAMP2 is then reexecuted and the nozzle contour is 
modified using the displacement thickness information contained 
on Unit 2. After the modified contour flowfield solution is 
completed, BLIMPJ is again executed to calculate the nozzle 
boundary layer using edge conditions taken at the wall of the 
modified nozzle flowfield solution. While in the BLIMPJ 
solution, a check is made once the boundary layer calculations 
have been completed. If the pressure at the lip of the nozzle 
varies by more than 10% from the pressure at the inside edge of 
the boundary layer in the inviscid solution (modified contour) 
then another pass is made through BLIMPJ with the edge conditions 
interpolated from the nozzle solution based on the location of 
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the boundary layer edge in the flowfield. These two passes are 
made within the last BLIMPJ execution and require no intervention 
by the user. This second pass is made to ensure that the 
boundary layer and inviscid flowfield properties match at the 
edge of the boundary layer. 

The MRE 15 engine specifications that were used as partial 
input for this analysis are contained in Table 1. Table 2 
presents the TRAN72 input file for generating the thermodynamic 
data input file for RAMP2 . Table 3 contains the input files for 
the RAMP2 code. The nozzle contour represented by a table of X, 

Y and flow angle values is contained in the RAMP2 input data 
shown in Table 3. The RAMP2 flowfield was calculated on the 
NASA/MSFC ED31 VAX. The VAX COM files for executing the codes 
are shown in Table 4 (TRAN72) and Table 5 (RAMP2) . 

The RAMP2 flowfield characteristics of the MRE 15 exhaust 
plume are shown in Figures 1-9. These figures present isovalue 
contours of Mach number, static pressure, temperature, pitot 
total pressures, density, velocity, mass flux, flow angle, and 
constant mass flow streamlines. The RAMP2 binary output file 
(Unit 3) is stored on the NASA/MSFC VAX and can be obtained by 
contacting the NASA COR for this study. 

3.2 MRE15 Source Flow Plume Model 

Recently, a model has been developed which fits the RAMP2 
flowfield results using a source flow approximation. This model 
called SFPGEN (Ref. 2) was developed to support a Space Station 
Contamination Model developed by SEA, Inc. for NASA (Ref. 6) . 

The source flow plume model generator is a code which 
determines the flow properties in the RAMP2 flowfield. The 
SFPGEN code is based on the radial lookup code (Ref. 7, 8). 
Modifications have been made to this code to determine the local 
properties in the already computed RAMP2 flowfield and then 
output a formatted file of a simple representation of the RAMP2 
flowfield. Properties are determined at constant distances from 
the exit plane as a function of angular distance from the plume 
axis. Ten of these source surfaces are determined at even 
increments of distance starting at the exit plane of the motor 
and terminating at the maximum distance on the axis in the RAMP2 
plume. RAMP2 should be executed out to about 80-100 exit radii 
where all the gas streamlines no longer change slope. At each of 
the 10 data surfaces the local/centerline ratios of density, 
total enthalpy and temperature are determined as a function of 
angle off the axis, out to the maximum expansion angle. The 
maximum expansion angle is passed to the source flow module from 
RAMP2 . In addition to the angular variation in properties, the 
code also generates an array which contains the flow property 
(density, velocity and total enthalpy) distribution along the 
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axis. Twenty five points are used at even axial increments, 
starting at the exit plane and terminating at the last source 
station. A single set of mass flow averaged species 
distributions is also determined at the exit plane of the nozzle. 
A formatted output file (Unit 16) is then generated which can be 
passed on as input data to the Space Station Contamination Model 
or PLIMP code. 

The source ^ flow plume model has the advantage over other 
source models, in that it has the capability to include total 
temperature effects in the exhaust plume due to oxidizer/fuel 
distribution and nozzle wall boundary layer effects. This model 
does assume a constant specific heat ratio and molecular weight. 
The specific heat ratio and species distribution which is used 
are mass flow averaged values at the exit plane. Engines which 
are typically used for orbital application have high area ratios 
where the exit plane specific heat ratio is very close to what 
the specific heat ratios would be in the plume. Additionally 
(neglecting 0/F ratio effects) the species at the exit plane are 
frozen and do not change in the plume. Previous comparisons of 
RAMP2 and SFPGEN show that SFPGEN results almost exactly match 
the RAMP2 calculations. 

The baseline SFPGEN code used 10 radial distributions and 25 
points on the axis to describe the flowfield. The radial 
distribution of flowfield properties is assumed constant from the 
nozzle exit plane up to the first source station. Therefore, in 
the near field (< 8 exit diameters) , the SFPGEN flowfield 
properties can be in error due to the large gradients which are 
present near the exit plane. 

Using this model during the checkout of PLIMP, the flowfield 
properties, as expected, in the near field differed enough to 
cause significant differences in predicted environments using the 
SFPGEN and RAMP2 plumes. To improve the near field results 
SFPGEN was modified to fit the plume with 25 source lines and up 
to 50 axis points. In addition, the source lines were made more 
concentrated in the near field through the use of a cosine 
distribution. The source lines also begin 1.5 exit radii from 
the exit plane. These modifications greatly increased the 
accuracy in the near field distribution while maintaining the 
fidelity in the far field. At high impingement angles, the 
boundary layer edge velocity (important for convective heating 
rate predictions) is very sensitive to local flow angle. The 
baseline SFPGEN code assumed that the local flow angle was the 
arctan (Y/X) . For the far field, this assumption is valid. 
However, for the near field (< 10 diameters) and beyond 80 
degrees where boundary layer effects are significant, this 
assumption can result in flow angles that are in error by as much 
as 7 degrees. To improve the flow angle prediction, flow angle 
was added to the local/centerline ratios at each source line. 
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This modification improved the flow angle prediction to + 1.5 
degrees. A detailed description of the SFPGEN program can be 
found in Reference 2. 

The major advantages of using the SFPGEN code to produce a 
flowfield are: the flowfield is portable between different 

computers; the flowfield requires considerably less computer 
storage than the RAMP2 flowfield; auxiliary (induced . 
environments) calculations are considerably faster using a source 
flow plume; the source flow plume can be extended to much larger 
distances than are practical for running the RAMP2 code; and the 
model can be easily used by other codes. 

The RAMP2 flowfield description is stored on a binary file 
which requires that the same computer type be used to generate 
the plume and perform application calculations. The SFPGEN 
flowfield output is a formatted ASCII file that can be easily 
transported from computer to computer without rerunning the RAMP2 
code. Thus, once a particular RAMP2 flowfield is generated, the 
SFPGEN flowfield can be transferred to other organizations very 
easily without the cost of rerunning the RAMP2 plume. 

A typical RAMP2 binary flowfield file requires on the order 
of 4-5 megabytes of mass storage. The SFPGEN flowfield requires 
about .06 megabytes of storage. Thus, the computer resource 
requirements are significantly reduced. 

Considerable time is spent by application codes (PLIMP, Ref. 
9, 10) in looking up flowfield properties in a RAMP2 flowfield. 
Timing tests performed using RAMP2 and SFPGEN flowfields show 
that PLIMP runs 3-4 times faster using a SFPGEN flowfield. For 
applications which require significant PLIMP calculations, 
significant reductions in computer resources can be realized. 

Reference 6 contains a description of the original routines 
which use this data file to determine local flow properties in 
the plume. These routines have been subsequently modified and 
copies of these routines, as well as the MRE 15 SFPGEN data file, 
are maintained on the NASA/MSFC VAX and can be obtained from the 
NASA - COR. These routines will be particularly useful to 
organizations that are interested in plume induced environments 
to payloads being deployed or retrieved by the OMV. 

Figures 10 - 16 show the flowfield characteristics of the 
MRE 15 SFPGEN plume. These figures present isovalve contours of 
Mach number, static pressure, temperature, pitot total pressure, 
density, velocity and mass flux, respectively. 

4.0 Conclusions 

This report provides a detailed description of the OMV 
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vehicle MRE 15 exhaust plume flowfield characteristics calculated 
with the RAMP2 and SFPGEN codes. Included in the report are 
brief descriptions of the two models, MRE 15 engine 
characteristics and RAMP2 input data files. The RAMP2 flowfield 
could be recalculated by other organizations using the 
information contained in this report. The MRE 15 SFPGEN 
flowfield can be readily used by other organizations who are 
interested in MRE 15 plume induced environments which require 
local flowfield properties which can be supplied using the SFPGEN 
MRE 15 model. 
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Table 1. MRE-15 Engine Specifications 


Inlet Pressure (psia) 300 

Nominal Thrust (lbf) 15 

Area Ratio 50 

Chamber Pressure (psia) 134.7 

Chamber Temperature (*F) 1700 

Throat Diameter (in) .296 ± .005 

Propellant Hydrazine 

Nozzle Contoured 

Nozzle Wall Temperature (Throat) (*F) 1700 

Nozzle Wall Temperature (Lip) (*F) 1300 
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Listing TRAN 7 2 VAX COM File For MRE15 Engine 


$ 

$ 

$ 

$ 

$ 

$ 


C03 3013.RAMPD 

AsIlGN MPr!? E15,DAT F0R005 
Assirw JS?iS2! 10 * DAT foroio 
JS fJ®! ^CTHERM.R2F FOR004 

PI S tm I ?LS EC,0UT FOR006 
RUN C033013 .RAMP3CECR2F 


Table 5 

Listing of RAMP 2 VAX COM File For MRE15 Engine 


v LU33013 .RAMP 3 

$ ASSIGN MRE15.0UT F0R006 
$ ASSIGN MRE15UN7.DAT FOR007 

$ ASSIGN MRE15UN8.DAT 

$ ASSIGN MRE15UN9.DAT 

ASSIGN MRE15UN10.DAT 

iCCTPlt kmni 


' , JJ/VJ 

$ ASSIGN MRE15UN3.DAT 

A r*nmr 


F0R008 

FOR009 

FOROIO 

FOR003 


r UK u U 3 

A?^r£ E S2P M,R2F MKE15UN4A.DAT 
ASSIGN MRE15UN4A.DAT F0R004 

A R R T MDn crnir * 


$ 

$ 

$ 


’ * ‘uviui u/il 

f £ff^ GN MRE1 5UN5A.DAT FOR005 
$ ^®F15UN1 .DAT F0R001 

ASSIGN MRE15UN2.DAT F0R002 
ASSIGN MRE15UN12.DAT F0R012 
t A ffJ GN MRE15UN13.DAT F0R013 
f MRE15 UN14.DAT F0R0I4 

$ ASSIGN MRE15UN17.DAT F0R017 
$ ASSIGN MRE15UN11.DAT F0R011 
$ RUN C033013 .RAMP3RAMP2P 
$ ASSIGN BLIMPJ.OUT FOR006 
$ RUN C033013 .RAMP3BLIMPJ 

t 5UN2 . DAT MRE15UN2A.DAT 

$ ASSIGN MRE15.0UT1 F0R006 

S £ffJ GN MRE1 5UN5B.DAT F0R005 
£ MRE15UN1A.DAT FOROOl 

$ RUN C033013.RAMP3RAMP2P 
£ HRE15UN3.DAT MRE15UN3B.DAT 
£ c^?t IGN BLIMp J*OUTl FOR006 
J C033013 .RAMP3BLIMPJ 

f fOPY HRE15UN1A.DAT MREUN1A . DAT1 
$ ASSIGN FOROOl . DAT FOROOl 
$ ASSIGN MRE15.0UT2 FOR006 
£ oS! 1 ? 11 MREl 5UN5C.DAT FOR005 
$ RUN C033013.RAMP3RAMP2P 
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Figure 1. Mach Number Contours for MRE15 Engine Exhaust Plume 
Calculated Using RAMP 2 
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Figure 2. Pressure Contours for MRE15 Engine Exhaust Plume 
Calculated Using RAMP2 
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Calculated Using RAMP2 
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Figure 4. Pitot Total Pressure Contours for MRE15 Engine 
Exhaust Plume Calculated Using RAMP 2 
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Figure 6. Velocity Contours for MRE15 Engine Exhaust Plume 
Calculated Using RAMP2 
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Figure 8. Flow Angle Contours for MREl 5 Engine Exhaust Plume 
Calculated Using RAMP2 
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Figure 9. Mass Flow Contours for MRE15 Engine Exhaust Plume 
Calculated Using RAMP2 
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Figure 10. Mach Number Contours for MRE15 Engine Exhaust Plume 
Calculated Using SFPGEN 
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Figure 11. Pressure Contours for MRE15 Engine Exhaust Plume 
Calculated Using SFPGEN 
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Figure 12. Temperature Contours for MRE15 Engine Exhaust Plume 
Calculated Using SFPGEN 
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Figure 13. Pitot Total Pressure Contours for MREl 5 Engine 
Exhaust Plume Calculated Using SFPGEN 
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Figure 15. Velocity Contours for MRE15 Engine Exhaust Plume 
Calculated Using SFPGEN 
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Figure 16. Mass Flux Contours for MRE15 Engine Exhaust Plume 
Calculated Using SFPGEN 
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1.0 Forward 

This document presents the results of work performed by 
SECA, Inc. for NASA, Marshall Space Flight Center under Contract 
NAS8-38243 . The NASA contracting officer representative for this 
study is Mr. P. R. Sulyma, ED33. 

2.0 Introduction 

This report presents the flowfield characteristics of the 
Orbital Maneuvering Vehicle (OMV) , GN2 thruster exhaust plume. 

The GN2 thruster is a 5 lbf cold gas (Nitrogen) engine which will 
be used as the secondary engine for the Reaction Control System 
(RCS) of the OMV. The GN2 thruster plume characterization is 
required in order to predict the GN2 thruster plume induced 
environments to the OMV itself, as well as to payloads it is 
retrieving or deploying. This report presents the GN2 thruster 
flowfield characteristics and describes the methodology which was 
used to characterize the plume. 

3.0 Discussion 

This section describes the GN2 thruster exhaust plume 
flowfield which was generated using the RAMP2 (Ref. 1) code and a 
new source flow model (Ref. 2). 

3.1 RAMP2 GN2 Thruster Flowfield 

The analysis of the GN2 nozzle and plume flowfields was 
performed using the Reacting and Multi-phase Computer Code (RAMP2 
Ref. 1) . RAMP2 is a method of characteristics based code used to 
solve inviscid supersonic flowfields. 

The purpose of the RAMP2 code is to solve a rocket exhaust 
flowfield in sufficient detail so that the plume may be used to 
adequately specify some design environment for systems which are 
influenced by the plume. For high altitude plumes, an adequate 
definition of the plume requires the inclusion of chemistry and 
nozzle boundary layer effects. The TRAN72 (Ref. 3) program and 
BLIMPJ (Ref. 4) codes are included in the RAMP2 code along with 
the RAMP2F flowfield code to treat chemistry and boundary layer 
effects . 

In order to use the original RAMP2 code (Ref. 5) to solve a 
high altitude plume, the following steps are required. First, 
the TRAN72 program input data is prepared and executed to 
generate a data tape describing the thermodynamic characteristics 
of the post-combustion gases. Next the RAMP2F flowfield data are 
prepared and the nozzle flowfield is solved using the TRAN72 
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program data tape as input. Then in order to adequately describe 
the nozzle boundary layer, the BLIMPJ code is executed using an 
input data file and flowfield tape generated by the RAMP2F nozzle 
solution. Finally, the exhaust plume is generated by using the 
nozzle solution and boundary layer solution to generate an exit 
plane start line that is used to initiate the plume solution. 

Thus, the generation of a high altitude plume could require up to 
four different executions of programs (TRAN72, RAMP2F, BLIMPJ, 
and RAMP2F) for the specification of the most detailed and 
accurate results. Physical input data are required only for the 
TRAN72 and first RAMP2F execution. All data required for the 
BLIMPJ code and second RAMP2F execution are generated internally 
to the program and/or communicated via temporary files. 

Recently RAMP2 has been modified (Ref. 1) which will better 
handle nozzles which operate at a low chamber pressure and/or 
have high area ratio nozzles. Motors of these types may have 
boundary layers which contain a significant percentage of the 
mass flow of the motor. Large boundary layers can interact with 
the inviscid portion of the flow and alter the flowfield from 
what would be calculated, assuming a pure inviscid calculation. 
The original version of RAMP2 which was mentioned above (Ref. 5) 
could not account for the interaction of the boundary layer and 
the inviscid flowfield. RAMP2 made a single pass through the 
nozzle and performed a single boundary layer solution using the 
BLIMPJ code. These solutions were then merged at the exit plane 
and plume solution was performed. 

The latest version of RAMP2 (Ref. 1) was used to calculate 
the GN2 thruster plume. This version allows a first order 
approximation of the interaction of the boundary layer and 
inviscid flowfield. The basic method involves making two passes 
through RAMP2 to calculate the nozzle flowfield and two passes 
through the BLIMPJ code to calculate the boundary layer. The 
first pass through the RAMP2 code calculates the nozzle flowfield 
using the original nozzle contour. After the first nozzle 
solution, the BLIMPJ code is executed using the conditions which 
exist at the nozzle wall (from the first nozzle solution) as the 
inviscid streamline boundary layer edge conditions. The results 
of this calculation including the boundary layer displacement 
thickness as a function of nozzle station are stored on FORTRAN 
Unit 2. RAMP2 is then reexecuted and the nozzle contour is 
modified using the displacement thickness information contained 
on Unit 2. After the modified contour flowfield solution is 
completed, BLIMPJ is again executed to calculate the nozzle 
boundary layer using edge conditions taken at the wall of the 
modified nozzle flowfield solution. While in the BLIMPJ 
solution, a check is made once the boundary layer calculations 
have been completed. If the pressure at the lip of the nozzle 
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varies by more than 10% from the pressure at the inside edge of 
the boundary layer in the inviscid solution (modified contour) 
then another pass is made through BLIMPJ with the edge conditions 
interpolated from the nozzle solution based on the location of 
the boundary layer edge in the flowfield. These two passes are 
made within the last BLIMPJ execution and require no intervention 
by the user. This second pass is made to ensure that the 
boundary layer and inviscid flowfield properties match at the 
edge of the boundary layer. 

While the latest version of the RAMP2 code was used to 
calculate the OMV GN2 thruster nozzle/plume flowfield, the 
multiple pass option was not utilized. The OMV GN2 nozzle has a 
low area ratio (10) nozzle that will result in negligable 
boundary layer effects on the inviscid flowfield. Additionally, 
it was not necessary to use the TRAN72 program to calculate the 
gas thermodynamics since the nitrogen gas can be treated as an 
ideal gas having a constant molecular weight (28) and specific 
heat ratio (1.4) 

The GN2 engine specifications that were used as partial 
input for this analysis are contained in Table 1. Table 2 
contains the input files for the RAMP2 code. The RAMP2 flowfield 
was calculated on the NASA/MSFC ED31 VAX. The VAX COM file for 
executing the codes is shown in Table 3 (RAMP2). 

The RAMP2 flowfield characteristics of the GN2 thruster 
exhaust plume are shown in Figures 1-9. These figures present 
isovalue contours of Mach number, static pressure, temperature, 
pitot total pressures, density, velocity, mass flux, flow angle, 
and constant mass flow streamlines. The RAMP2 binary output file 
(Unit 3) is stored on the NASA/MSFC VAX and can be obtained by 
contacting the NASA COR for this study. 

3.2 GN2 Thruster Source Flow Plume Model 

Recently, a model has been developed which fits the RAMP2 
flowfield results using a source flow approximation. This model 
called SFPGEN (Ref. 2) was developed to support a Space Station 
Contamination Model developed by SEA, Inc. for NASA (Ref. 6). 

The source flow plume model generator is a code which 
determines the flow properties in the RAMP2 flowfield. The 
SFPGEN code is based on the radial lookup code (Ref. 7, 8). 
Modifications have been made to this code to determine the local 
properties in the already computed RAMP2 flowfield and then 
output a formatted file of a simple representation of the RAMP2 
flowfield. Properties are determined at constant distances from 
the exit plane as a function of angular distance from the plume 
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axis. Ten of these source surfaces are determined at even 
increments of distance starting at the exit plane of the motor 
and terminating at the maximum distance on the axis in the RAMP2 
plume. RAMP2 should be executed out to about 80-100 exit radii 
where all the gas streamlines no longer change slope. At each of 
the 10 data surfaces the local/centerline ratios of density, 
total enthalpy and temperature are determined as a function of 
angle off the axis, out to the maximum expansion angle. The 
maximum expansion angle is passed to the source flow module from 
RAMP2 . In addition to the angular variation in properties, the 
code also generates an array which contains the flow property 
(density, velocity and total enthalpy) distribution along the 
axis. Twenty five points are used at even axial increments, 
starting at the exit plane and terminating at the last source 
station. A single set of mass flow averaged species 
distributions is also determined at the exit plane of the nozzle. 
A formatted output file (Unit 16) is then generated which can be 
passed on as input data to the Space Station Contamination Model 
or PLIMP code. 

The source flow plume model has the advantage over other 
source models, in that it has the capability to include total 
temperature effects in the exhaust plume due to oxidizer/fuel 
distribution and nozzle wall boundary layer effects. This model 
does assume a constant specific heat ratio and molecular weight. 
The specific heat ratio and species distribution which is used 
are mass flow averaged values at the exit plane. Engines which 
are typically used for orbital application have high area ratios 
where the exit plane specific heat ratio is very close to what 
the specific heat ratios would be in the plume. Additionally 
(neglecting 0/F ratio effects) the species at the exit plane are 
frozen and do not change in the plume. Previous comparisons of 
RAMP2 and SFPGEN show that SFPGEN results almost exactly match 
the RAMP2 calculations . 

The baseline SFPGEN code used 10 radial distributions and 25 
points on the axis to describe the flowfield. The radial 
distribution of flowfield properties is assumed constant from the 
nozzle exit plane up to the first source station. Therefore, in 
the near field (< 8 exit diameters) , the SFPGEN flowfield 
properties can be in error due to the large gradients which are 
present near the exit plane. 

Using this model during the checkout of PLIMP, the flowfield 
properties, as expected, in the near field differed enough to 
cause significant differences in predicted environments using the 
SFPGEN and RAMP2 plumes. To improve the near field results 
SFPGEN was modified to fit the plume with 25 source lines and up 
to 50 axis points. In addition, the source lines were made more 
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concentrated in the near field through the use of a “sine 

accuracy in the near field impingement angles, the 

boundary layei^edge ^ity^Crtan? J-onvective^eating 
rate predictions) is very local flow angle was the 

baseline SFPGEN code assumed that «“ 1“ tion J valid. 

However for' the°near e field « li diameters) and beyon^eo 

degrees where boundary layer effects s ign £ b much 

T^improve^the floi angle 

was added to the local/centerline ratios at source 1 ^ ^ 

degrees^ '^detailed^description^ of III SF?gL program can be 
found in Reference 2. 

The major advantages of using the SFPGEN to Produce a 

the flowfield is portable between different 

model can be easily used by other codes. 

Th . -RAMP2 flowfield description is stored on a binary file 

which requires that the same computer be use ^ h g° s ^pQEJ a 6 

the plume and perform application cacu ®tion ^ be easily 

1? pSen f lowf iel^can'be transferred to other organisations very 
easily without the cost of rerunning the RAMP2 plume. 

A typical ^^^-ty o flo„field 

about 5 . 06 9 megabjtes of storage. Thus, the computer resource 
requirements are significantly reduced. 

considerable time is spent by application codes f ^ 

dt tmp runs 3-4 tunes faster usin^ a . 


Reference 6 


contains a description of the original routines 
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which use this data file to determine local flow properties in 
the plume. These routines have been subsequently modified and 
copies of these routines, as well as the MRE 15 SFPGEN data file, 
are maintained on the NAS A/MS FC VAX and can be obtained from the 
NASA - COR. These routines will be particularly useful to 
organizations that are interested in plume induced environments 
to payloads being deployed or retrieved by the OMV. 

Figures 10 - 16 show the flowfield characteristics of the 
GN2 SFPGEN plume. These figures present isovalve contours of 
Mach number, static pressure, temperature, pitot total pressure, 
density, velocity and mass flux, respectively. 

4.0 Conclusions 

This report provides a detailed description of the Oribital 
Maneuvering Vehicle GN2 thruster exhaust plume flowfield 
characteristics calculated with the RAMP2 and SFPGEN codes. 
Included in the report are brief descriptions of the two models, 
GN2 thruster characteristics and RAMP2 input data files. The 
RAMP2 flowfield could be recalculated by other organizations 
using the information contained in this report. The GN2 
flowfield can be readily used by other organizations who are 
interested in GN2 plume induced environments which require local 
flowfield properties which can be supplied using the SFPGEN GN2 
model . 
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Table 1. GN2 Thruster Specifications 


Inlet Pressure (psia) 

276-315 

Nominal Thrust (lbf) 

5 

Area Ratio 

10 

Chamber Pressure (psia) 

63 

Chamber Temperature (°F) 

530 

Throat Diameter (in) 

.25 ± .001 

Propellant 

Nitrogen 

Nozzle 

15° Conical 

Adiabatic Nozzle Wall 

Nozzle Length (Throat to Exit) (In) 

1.0241 
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Figure 1. Mach Number Contours for GN2 Thruster Exhaust Plume 
Calculated Using RAMP2 
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Figure 2. Pressure Contours for GN2 Thruster Exhaust Plume 
Calculated Using RAMP2 
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iqure 3. TemperaLure Contoum 1 
Calculated Usincj RAMP 2 
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Calculated Using RAMP2 
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Figure 5. Density Contours for GN2 Thruster Exhaust Plume 
Calculated Using RAMP 2 
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Calculated Using RAMP2 
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Figure 7. Mass Flux Contours for GN2 Thruster Exhaust Plume 
Calculated Using RAMP 2 
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l-’itjuie 8. Flow Antj I *■' font fMirs f or CM/! Thruster Exhaust 
Calculated Usinq RAMP2 
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i'xcjurc 9. Mass Flow Contours lor GN2 Thruster Fxhaust Plume 
Calculated Usinq RAMP2 
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Mach Number Contours lor CM2 Thruster Exhaust Plume 
Calculated Usinq SFPGEN 
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I’rossure Contours lor (>N 2 I'hv uslor Exhaust- I 1 unu 
Calculated Using SFPGEN 
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juro 1 / . Temperat ure' Ontoui’s f'>r ON? Thruster Exhaust. N umr 
Calculated Using SFPGEN 



□MU COLD GRS NOZZLE/PLUME GNP PC=63 


SECA-TR-9 0-1 3 




<S> 


ru 

ru 

m 

m 

XT 


XT 

LH 

LD 

LO 

CD 


<5> 

<s> 

<s> 

s> 

<s> 

<s> 

<£> 

<5> 

<s> 

<5> 

<S> 

<S> 

<s« 

<S> 


4* 

4- 

1 

i 

i 

i 

i 

i 

1 

t 

i 

i 

1 

i 

CE 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

L L- 

<S> 

<s> 

<S> 

<S> 

S> 

<S> 

<S> 

<5> 


<S> 

<S> 


<S> 

<S> 

If) 

<s> 

<s> 

<5> 

<S> 

S> 

<S> 


<S> 

<s> 

<S> 

<S> 

<S> 


<S> 

CL 


ru 

<s> 

fH 


LD 

XT 

00 

m 

XT 

LD 

LD 

cr> 


w 

Ln 

V 

ru 

CD 

r- 

CD 

<S> 

<S> 

xr 

LO 

ru 

LO 

r- 

ru 

CL 


OJ 

m 

CD 


LD 

ru 

CO 

m 


r- 

CO 


cn 

i 

r n 

<5> 

s> 


<S> 

<s> 

<S> 

<s> 

<s> 

<s> 

<S> 

<s> 

<s> 

«s> 

<s> 

<E 

CD 

cc 

oo 

u 

o 

LU 

Ll 

LD 

X 

»— i 



— 1 





^ LO ® LD ® LO ® 

ro ru ru ^ 

-u 




JZ 

E- 


23 


Pitot To t . 7 i ] Pi'-ssuro Contours for GN 2 

Exhaust Plume Calculated Using SFPGEN 
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Calculated Usiruj Si PGEN 
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Velocity Contours lor CM2 Thrust or Exhaust Plurru 
Calculated Using SEPGEN 



